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Hidden topological angles and Lefschetz thimbles 
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We demonstrate the existence of hidden topological angles (HTAs) in a large class of quantum field theories 
and quantum mechanical systems. HTAs are distinct from theta-parameters in the lagrangian. They arise as 
invariant angle associated with saddle points of the complexified path integral and their descent manifolds 
(Lefschetz thimbles). Physical effects of HTAs become most transparent upon analytic continuation in ny to 
non-integer number of flavors, reducing in the integer ny limit to a Z 2 valued phase difference between dominant 
saddles. In = 1 super Yang-Mills theory we demonstrate the microscopic mechanism for the vanishing of 
the gluon condensate. The same effect leads to an anomalously small condensate in a QCD-like SU ( N ) gauge 
theory with fermions in the two-index representation. The basic phenomenon is that, contrary to folklore, the 
gluon condensate can receive both positive and negative contributions in a semi-classical expansion. In quantum 
mechanics, a HTA leads to a difference in semi-classical expansion of integer and half-integer spin particles. 


Introduction. Providing a non-perturbative continuum defi¬ 
nition of the path integral in quantum field theory is a chal¬ 
lenging but important problem [1]. There is growing ev¬ 
idence that, if an ordinary integral or a path integral ad¬ 
mits a Lefschetz-thimble decomposition [2, 3] or resurgent 
transseries expansion [4-12] then either of these methods 
gives this long-sought non-perturbative definition. If this is 
indeed the case, then we expect that these new methods will 
provide new and deep insight into quantum field theory and 
quantum mechanics formulated in terms of path integrals. In 
this article we introduce a new phenomenon of this kind, the 
appearance of hidden topological angles (HTAs). 

The main prescription associated with the Lefschetz- 
thimble decomposition or the resurgent expansion is the fol¬ 
lowing: Even if an ordinary integral or a path integral is 
formulated over real fields, the natural space that the criti¬ 
cal points (saddles) p a live in is the complexification of the 
original space of fields. However, the dimension of the crit¬ 
ical point cycles J7 a is that of the original space, or half that 
of the complexified field space. For example, for an ordi¬ 
nary integral over A-dimensional real space, this procedure is 
R n — C N —» L n , where = £ 0 ft 0 J7o and dim^Jfc) = N. 
For N = 1, this is the well-known steepest descent (stationary 
phase) approximation. 

To each critical point p c of the complexified action one 
attributes an action, with real and imaginary parts, and with 
“weight” e~ Sa . The imaginary part of the action, Im5 0 is 
an invariant angle associated with the critical point p c and 
its descent manifold J7 0 . If there are critical points with the 
identical real part of the action Re5 0 , but different imaginary 
parts Im5 C , then there may be subtle effects. Indeed, Wit¬ 
ten recently studied in [3] the analytic continuation of Chern- 
Simons theory to non-integer values of the coupling k , find¬ 
ing subtle cancellations among dominant saddle field config¬ 
urations in the integer k limit, so that the sub-dominant sad¬ 
dle gives the main physical contribution. In this work, we 
show that the effect observed by Witten is not an exotic phe¬ 
nomenon, but that it is possibly quite ubiquitous, and that 
it is responsible for a variety of interesting physical effects 
in quantum field theories and quantum mechanics, in which 
the analytic continuation is now to non-integer “coupling” ny. 


which is the number of fermionic flavors for integer values. 
We also show that the effect is more non-trivial than a simple 
cancellation between dominant saddles. Indeed, the effect de¬ 
pends on the observable, and the dominant saddles may can¬ 
cels in certain observables, but contribute to others. 

In a field theory with a topological 0-angle in the La¬ 
grangian, subtle effects may arise at certain values of the 0- 
angle [13-16]. For example, at 0 = n in 51/(2) gauge the¬ 
ory there is a cancellation of leading order saddle contribu¬ 
tion to the mass gap [14]. In this work we study a more ex¬ 
otic phenomenon, which is due to a hidden topological angle 
not explicitly present in the lagrangian. We define a HTA as 
the phase associated with a saddle point in the complexified 
field space. Below, we will provide examples of this phe¬ 
nomenon in = 1 super Yang-Mills theory (SYM), certain 
QCD-like theories, and the quantum mechanics of a particle 
with spin. We also note that a HTA is different from the dis¬ 
crete theta angles discussed recently [17], which comes about 
as one changes the global gauge group. In contrast, HTAs are 
present for any gauge group. In the examples discussed below 
we find that for integer values of the number of fermions ny 
there is a Z 2 hidden topological structure. 

A prototype in ordinary integration. An elementary exam¬ 
ple that provides some intuition for field theory is the follow¬ 
ing. Consider the analytic continuation of the Bessel func¬ 
tion to non-integer order, and describe the contour that ap¬ 
pears in the integral representation in terms of Lefschetz thim¬ 
bles. In one complex dimension the Lefschetz thimble is de¬ 
fined as a stationary phase manifold: Im[5(w) — S(w n )] = 0 
where w n is a critical point on the contour. The integral is 
I(k,X) = f Cw dwe 2 ^ smh ^ +kw for complex k , X. In a certain 
regime of the analytic continuation, discussed in [3], the inte¬ 
gral can be expressed in terms of three cycles, %,i= 1,2,3 as¬ 
sociated with saddles p ? , so that C w = J\ + J 2 + _? 3 , see Fig. 1. 
The sum of the three thimble contribution gives 

^ l+e 2ni{k+\)y-S l+e -S 2 (1) 

where \e~ Sx | = \e~ Si | A> e~ Sl , i.e, pi and p 3 are dominant over 
p 2 . However, they have a relative phase, and the contribution 
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of these two dominant saddles cancel each other exactly for 
integer k. The mechanism described above is an intuitive ex¬ 
ample of a mechanism operative in Chern-Simons theory by 
using analytic continuation, providing confidence for the util¬ 
ity of the idea of analytic continuation of path integrals. Other 
examples are discussed in [18-20]. We will perform a similar 
analytic continuation in ttf, the number of fermion flavors in 
the theory. 



FIG. 1: The blue areas show “good regions” in which the integrand 
falls sufficiently rapidly at infinity to guarantee convergence. The red 
dots give the locations of the saddle points, and the blue contours are 
the Lefschetz thimbles. If the boundary of integration is (—°°, — 00 + 
271/), then the Lefschetz decomposition is J\ + J 2 + J 6 - Here, pi 
and P 3 are equally dominant saddles over P 2 , but there is an over-all 
phase difference between the dominant saddles leading to a subtle 
cancellation for integer k. 


Picard-Lefschetz equation and invariant angles. The 

definition of the Lefschetz thimble based on stationary 
phase, ImfS^w) — S(w n )] = 0, is only satisfactory for a one¬ 
dimensional integral (it provides one real condition on a one- 
complex dimensional space). In n complex dimensions, where 
n = 00 corresponds to field theory or quantum mechanics, this 
condition defines a co-dimension one (real dimension 2n — 1 
space), which is not the desired n real dimensional space. 
Instead, one needs n real conditions to define the thimble. 
Guided by these observations, Witten used complex gradient 
flow equations, the Picard-Lefschetz equations, to describe 
the Lefschetz-thimbles. In a theory with a field cp and action 
S( (p), this amounts to 

3cp _ 8 S dtp 8 S 

dx 8cp ’ dx Sep ’ 

where x is the flow time. Using (2) and the chain rule, 

aim [S] 1 fdsdcp 8Sdcp\ 

dx 2/ V Sep 0x Sep dx ) ’ 

meaning that Im[S(<|>)] = ImfS^)] is invariant under the flow. 
In a quantum field theory (QFT), or in quantum mechanics 


(QM), in which semi-classical saddle proliferates (an exam¬ 
ple is the instanton gas), Im[£((!)„)] will appear as a genuinely 
new phase in the effective field theory. This is the HTA phe¬ 
nomenon. 

The integration in the complexified field space is infinite 
dimensional. In the background of non-perturbative saddles, 
this space usually factorizes into finite dimensional zero and 
quasi-zero modes directions and infinite dimensional gaussian 
modes. The HTA can be calculated by an exact integration 
over the complexified finite dimensional quasi-zero mode di¬ 
rections in the field space, dictated by the finite dimensional 
version of the Picard-Lefschetz theory. 

Hidden topological angle in 4d = 1 SYM: Consider 
C\C = 1 SYM on M 3 x S[, where S[ is a circle with period 
L. We use supersymmetry preserving boundary conditions 
and take the small L limit in order to be able to use semi- 
classical methods. According to the trace anomaly relation 
the gluon condensate (^tr F 2 V ) determines the vacuum energy: 

“Evac = (£2|7oo|£2) = 5 (£212^1 £2) = i^|) (trF^}. This implies 
that the gluon condensate can serve as an order parameter for 
supersymmetry breaking. The vacuum energy density, and 
hence the condensate, vanishes to all orders in perturbation 
theory in supersymmetric theories. Since supersymmetry is 
known to be unbroken, the gluon condensate must be zero 
non-perturbatively as well. In the semi-classical limit this re¬ 
sult appears mysterious, because all contributions appear to 
be positive. The reason is that in euclidean space the fermion 
determinant is positive definite, and trF 2 v is also positive def¬ 
inite. This implies that the gluon condensate is the average of 
a positive observable with respect to a positive measure [21]. 
Then, how does the vanishing of the (tr F 2 V ) take place from a 
semi-classical point of view? 

We address this question in the regime of small, but finite 
radii on M 3 x S[. To do so, recall the Euclidean realization 
of the vacuum of the theory on small M 3 x S[, depicted in 
Fig. 2 for the center-symmetric point of the Wilson line on 
the Coulomb branch. The vacuum is, primarily, a dilute gas 
of semi-classical one- and two-events: monopole-instantons 
[22-25] and bions [5, 26-28]. These are: 

a) monopole-instantons, Mi = e~ s °(aj-X) 2 , 

b) magnetic bions, (2 \j = [MjM j\ = e~ 2S ° ..., 

c) neutral bions, (Bu = [Mi Mi] = e ~ 2S ° +m .... 

where oq, / = 1 ,...,N are simple roots complemented with the 
affine root a#, and ^ and %i are non-vanishing VA^- < 0, and 
MAu > 0 entries of the extended Cartan matrix, respectively. 
The monopole action is So = • For small L the coupling 

is small, the action is large, and fluctuations are suppressed. 
The IN fermion zero modes of the 4d instanton are distributed 
uniformly as (2,2,..., 2) to monopoles Mi. 

At leading order 0(e ~ s °) in the semi-classical limit, each 
monopole-instanton has two fermion zero modes and there- 
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FIG. 2: A snap-shot of the euclidean vacuum of = 1 SYM on 
small R 3 x Both neutral and magnetic bions carry action 2 Sq, but 
their contribution to gluon condensate cancels exactly because of the 
presence of a HTA, a 7t-phase difference between the two saddles. 


fore they do not contribute to the gluon condensate. Two- 
defects do contribute to the gluon condensate. For the sake of 
making the analogy with the toy example (1) explicit, let us 
consider analytic continuation away from nf = 1. The den¬ 
sity of both types of 2-defects is the same, of order 0(e~ 2S °). 
However, there is an extra (4 nf — 3)n phase (invariant angle) 
associated with the neutral bion saddle/thimble: 

Arg(j?®..) = Arg(j/® iy ) + (4 n f -3)n. (4) 

Consequently, in contrast to the folklore regarding the positiv¬ 
ity of the gluon condensate, the contributions of the two types 
of 2-defects to the gluon condensate cancel: 

^ 4 (^ tr ^uv) = 0 x n Mi + + e ^ Anf = 0• (5) 

at a physical integer value of the parameter, nf = 1, similar to 
(1). This is the microscopic mechanism for the vanishing of 
the gluon condensate as well as the vacuum energy in 7\£ = 1 
SYM. The two contributing bion-thimbles are charged oppo¬ 
sitely under the Zi ? TS , and cancel each other out. 

The difference with respect to the toy integral and the can¬ 
cellation in analytically continued Chern-Simons theory is the 
fact that this cancellation is observable dependent. In fact, the 
combination of the neutral and magnetic bions, despite giving 
vanishing contribution to gluon condensate, is responsible for 
the formation of a mass gap. To see this consider the effective 
lagrangian for the low energy bosonic modes. As an example, 
we will use SU (2) gauge theory. Let denote the fluctuation 
of the Wilson line around the center symmetric minimum and 
o denote the dual photon. The bosonic potential induced by 


2-defects is [29] 

V(CJ,C|)) = -(#12+ #21 +#11 + #22) 

~ e~ 2S ° (— cos2g — e m cosh 24 >) . (6) 

We observe that the factor e l7Z responsible for the vanish¬ 
ing (jjtrF 2 ,) is also responsible for the (positive and un¬ 
suppressed) mass gap of the (^-fluctuations and stabilizes 
center-symmetry. The HTA explains both the vanishing of the 
gluon condensate and the non-tachyonic nature of fluctuations 
of the Polyakov line. It is also not particular to supersymmet¬ 
ric theory, as we discuss next. 

QCD(AS/S): In a typical confining asymptotically free 
SU(N) gauge theory, the “natural” scaling of the (properly 
normalized) gluon condensate is 0(N°): (^tr F 2 V ) oc N° A 4 . It 
is natural to expect that the vanishing of the gluon condensate 
is special to the supersymmetric theory. This is not the case. 
There exists an exact large-TV orientifold/orbifold equivalence 
between = 1 SYM and QCD(AS/S) [30], proven in [31]. 
Here, AS/S refers to fermions in anti-symmetric/symmetric 
two-index representations. The largeW equivalence implies 
that the gluon condensate in QCD(AS/S) is zero at leading 
order in the 1 /N expansion, and must scale as [32]: 

<>+> QCD(AS/S) = +f (7) 

in sharp contrast with the “natural” value. This result is 
counter-intuitive, but it is a rigorous consequence of the large 
N equivalence. However, as in the supersymmetric case, there 
is no known semi-classical explanation. 

Again, we can understand the result based on the presence 
of HTAs. To achieve this, we use the framework of deformed 
Yang-Mills theory, and add AS representation fermions (a 
similar analysis holds for QCD(S)). In QCD(AS), the Atiyah- 
Singer index theorem implies that the number of fermion zero 
modes of a 4d instanton is 2N — 4. There are N types of 
monopole-instantons, with the number of fermion zero modes 
distributed as (2,2,... ,2,0,0) in a center-symmetric back¬ 
ground. The difference with respect to = 1 SYM is that 
2 out of N monopole-instantons do not possess fermi zero 
modes. Therefore, at leading order, 0(e ~ s °), in the semi- 
classical expansion, N — 2 monopoles do not contribute to the 
gluon condensate and only two do, giving a positive contribu¬ 
tion proportional to 1 /N. At second order in the semiclassical 
expansion, 0(e ~ 2S °), there are magnetic and neutral bions that 
can contribute to gluon condensate. Their contribution cancels 
at leading order in N , analogous to SYM, leading to (7). 

In the older literature on QCD [21, 33, 34], it was assumed 
that in the semi-classical limit gluon condensate is propor¬ 
tional to the instantons density. This was based on the ra¬ 
tionale that a single instanton contributes a finite and positive 
amount, ^2 / tr , and that the condensate can be at¬ 

tributed to 4d instantons with a positive weight, (^tr F 2 V ) oc 
nj. In the calculable small S[ regime, we see that this is in¬ 
correct in at least two ways: i) Instantons are sub-leading, i.e. 
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0(e~ N ), ii) The weight of the saddles can be both positive (de¬ 
creasing energy) and negative (increasing energy). Our work 
is the first example in which a contribution to condensate has 
a negative component. 

Quantum mechanics: In order to show the generality of hid¬ 
den topological angles, we also consider the quantum me¬ 
chanics of a particle with position x(t) and internal spin (\) nf . 
The euclidean Lagrangian is that of a bosonic field x(t) cou¬ 
pled to rif fermionic fields \|/*: 

Le=(\x 2 + \ (W') 2 + i(tjWi + rwO) (8) 

For rif = 1, this theory is supersymmetric [35]. If we choose 
W (x) to be a periodic function, for example W (x) = cosx, we 
may identify v = x + 2n as the same physical point, corre¬ 
sponding to a particle on a circle, rather than in an infinite 
lattice. The system contains two types of instantons, 

I\ : [0 —»• 7t], h : [0 -s- -rc]. (9) 

Here, h is an instanton (not an anti-instanton), because it sat¬ 
isfies the same BPS equation (or gradient flow equation, if W 
is viewed as a Morse function) as I \. 

Because of spin, instantons do not contribute to the vacuum 
energy. A non-vanishing contribution arises from correlated 
two-events. This parallels the 4d field theory on M? x S' 1 , and 
provides a simple system in which the effect of the HTA is 
not contaminated by first order instanton effects. Following 
[36, 37], we find that the amplitudes of the two-events are 
given by [I\Ii] = [I 2 I 2 ] e inn fe ~ 2Sl , and [/iT 2 ] = [J 2 / 1 ] oc 
e ~ 2Sl . Due to the difference in the invariant angles between 
the two-saddles/thimbles we find 

Ar §U[/Jj]) = Ar g(^[/,T 2 ]) + n f K ■ CIO) 

The non-perturbative contribution to the ground state energy 
takes the form: 

A£” p = (-2 - 2 e im f)e~ 2s ' (11) 

While the [/ 1 / 2 ] molecule behaves in the expected manner 
and decreases the ground state energy, the [I\ I\] molecule is 
sensitive to the HTA governed by the spin and increases the 
ground state energy for odd rif (half-integer spin) while de¬ 
creasing it for even-fty (integer spin). In the case rif = 1, (10) 
is the microscopic reason for the non-perturbative vanishing 
of the ground state energy. We note that this system, despite 
having Witten index zero [38], % = tr[(—1) F ] = 0, has unbro¬ 
ken supersymmetry, and two supersymmetric ground states. 

There are two additional interesting features of this system. 
The first is related to the fact that one can introduce a topolog¬ 
ical angle into the Lagrangian, / x Jx, and unlike the case 
of supersymmetric gauge theory, the 0 -angle is physical, and 
alters the spectrum of the theory. Since I w =0, the supersym¬ 
metry of this system is fragile. The vacuum energy can be 
written as 

AZsqP = (— 2 cos© — 2 e im f)e~ 2Sl ( 12 ) 


which, for the supersymmetric theory (rif = 1 ), takes the 
form: 

A£q P = (—2cos0 + 2)e -2iS/ = 4 sin 2 ®e~ 2Sl > 0 (13) 

meaning supersymmetry is dynamically broken for 0 / 0 . 
Note that the energy remains positive semi-definite, which is 
a consequence of the supersymmetry of the Hamiltonian. The 
physical reason for E > 0 in the case 0 ^ 0 is that the 0- 
angle is equivalent to feeding momentum into the system. Be¬ 
cause of supersymmetry, bosonic/fermionic ground state pair 
is lifted simultaneously by the insertion of momentum, lead¬ 
ing to a non-vanishing ground state energy. 

The second unusual feature is that the theory has Witten in¬ 
dex Iw = 0 for any value of 0 , but that the reason for Iw = 0 
differs in the two cases. For 0 = 0, we get % = 1 — 1 = 0, 
where the two contributions arise from the bosonic/fermionic 
sectors of the Hilbert space, and supersymmetry is unbroken. 
In the second case, 0 / 0, we get Iw = 0 — 0 = 0 and super- 
symmetry is broken. 

One may speculate that the invariant angles are related 
to Berry phases [39], realized in terms of Euclidean sad¬ 
dles, at least in the case of quantum mechanics. Since 
(\) nf = 0^mult(5')5', where mult(S') is the multiplicity, we 
can rewrite the path integral over the Grassmann variables as 
spin path integrals [40] Z = J^-mult^Z® where 

Z® = jDxD(cos6)D§ e ~ LE+i ^ IMx (14) 

L e = ±(x 2 + (W') 2 )+SW"cose + iS(l-cosG)ty 

where (0,4>) G S 2 parameterize the Bloch sphere. There are 
two spin dependent interactions, a “magnetic field” W"-spin 
coupling, and the Wess-Zumino term, or Berry phase action. 
In this language (10) should be replaced by ArgC?^ 7 ^) = 
Arg(j/j 7 i j 2 j) +2Sn. This distinguishes half-integer and inte¬ 
ger spin particles, similar to anti-ferromagnets in one-spatial 
dimension [41], leading to qualitative differences. 

Conclusion: We have provided several examples of Z 2 hid¬ 
den topological angles, associated with the saddle point man¬ 
ifolds that appear in the complexified path integral, and have 
shown that these angles lead to crucial physical effects. We 
anticipate that HTAs will have crucial impact on the semi- 
classical analysis of many interesting quantum field theories 
and quantum mechanical systems. Our examples also show 
that in an attempt to perform lattice simulations using Lef- 
schetz thimbles, e.g., [42-45], all thimbles whose multipliers 
are non-zero must be carefully summed over to correctly cap¬ 
ture the dynamics of the theory. 

Acknowledgments. We thank G. Dunne, G. ’t Hooft, P. 
Argyres, and E. Witten for useful comments and discus¬ 
sions. We acknowledge support from DOE grants DE-FG02- 
03ER41260 and DE-SC0013036. 


5 


* Electronic address: abehtas@ncsu.edu 
^ Electronic address: tsulejm@ncsu.edu 
^ Electronic address: tmschaef@ncsu.edu 
§ Electronic address: unsal.mithat@gmail.com 

[1] G. ’tHooft, Subnucl.Ser. 15, 943 (1979). 

[2] E. Witten, (2010), arXiv: 1009.6032 [hep-th] . 

[3] E. Witten, (2010), arXiv: 1001.2933 [hep-th]. 

[4] G. Y. Dunne and M. Unsal, Phys.Rev. D87, 025015 (2013), 
arXiv: 1210.3646 [hep-th]. 

[5] P. C. Argyres and M. Unsal, JHEP 1208, 063 (2012), 
arXiv: 1206.1890 [hep-th]. 

[6] G. Y. Dunne and M. Unsal, JHEP 1211, 170 (2012), 
arXiv: 1210.2423 [hep-th]. 

[7] G. V. Dunne and M. Unsal, Phys.Rev. D89, 041701 (2014), 
arXiv: 1306.4405 [hep-th]. 

[8] A. Cherman, D. Dorigoni, and M. Unsal, (2014), 
arXiv: 1403.1277 [hep-th]. 

[9] G. V. Dunne and M. Unsal, Phys.Rev. D89, 105009 (2014), 
arXiv: 1401.5202 [hep-th]. 

[10] I. Aniceto and R. Schiappa, (2013), arXiv: 1308.1115 [hep-th]. 

[11] T. Misumi, M. Nitta, and N. Sakai, JHEP 1406, 164 (2014), 
arXiv: 1404.7225 [hep-th]. 

[12] T. Misumi, M. Nitta, and N. Sakai, (2014), arXiv: 1409.3444 
[hep-th]. 

[13] E. Witten, Nucl.Phys. B507, 658 (1997), arXiv:hep-th/9706109 
[hep-th] . 

[14] M. Unsal, Phys.Rev. D86, 105012 (2012), arXiv: 1201.6426 
[hep-th]. 

[15] A. Bhoonah, E. Thomas, and A. R. Zhitnitsky, Nucl.Phys. 
B890, 30 (2014), arXiv: 1407.5121 [hep-ph]. 

[16] M. M. Anber, Phys.Rev. D88, 085003 (2013), arXiv: 1302.2641 
[hep-th]. 

[17] D. Gaiotto, A. Kapustin, N. Seiberg, and B. Willett, (2014), 
arXiv: 1412.5148 [hep-th]. 

[18] T. Kanazawa and Y. Tanizaki, (2014), arXiv: 1412.2802 [hep- 
th]. 

[19] Y. Tanizaki, (2014), arXiv: 1412.1891 [hep-th]. 

[20] Y. Tanizaki and T. Koike, Annals Phys. 351, 250 (2014), 
arXiv: 1406.2386 [math-ph]. 

[21] M. A. Shifman, A. Vainshtein, and V. I. Zakharov, Nucl.Phys. 
B147, 385 (1979). 


[22] K.-M. Lee and P. Yi, Phys.Rev. D56, 3711 (1997), arXiv:hep- 
th/9702107 [hep-th] . 

[23] K.-M. Lee and C.-h. Lu, Phys.Rev. D58, 025011 (1998), 
arXiv:hep-th/9802108 [hep-th] . 

[24] T. C. Kraan and P. van Baal, Nucl.Phys. B533, 627 (1998), 
arXiv:hep-th/9805168 [hep-th] . 

[25] N. M. Davies, T. J. Hollowood, and Y Y Khoze, J.Math.Phys. 
44, 3640 (2003), arXiv:hep-th/0006011 [hep-th] . 

[26] M. Unsal, Phys.Rev. D80, 065001 (2009), arXiv:0709.3269 
[hep-th]. 

[27] E. Poppitz and M. Unsal, JHEP 1107, 082 (2011), 
arXiv: 1105.3969 [hep-th]. 

[28] M. M. Anber and E. Poppitz, JHEP 1106, 136 (2011), 
arXiv: 1105.0940 [hep-th]. 

[29] E. Poppitz, T. Schafer, and M. Unsal, JHEP 1210, 115 (2012), 
arXiv: 1205.0290 [hep-th]. 

[30] A. Armoni, M. Shifman, and G. Veneziano, Nucl.Phys. B667, 
170 (2003), arXiv:hep-th/0302163 [hep-th] . 

[31] P. Kovtun, M. Unsal, and L. G. Yaffe, Phys.Rev. D72, 105006 
(2005), arXiv:hep-th/0505075 [hep-th] . 

[32] M. Unsal, Phys.Rev. D76, 025015 (2007), arXiv:hep- 
th/0703025 [HEP-TH] . 

[33] C. G. Callan, R. F. Dashen, and D. J. Gross, Phys.Rev. D17, 
2717 (1978). 

[34] T. Schafer and E. V. Shuryak, Rev.Mod.Phys. 70, 323 (1998), 
arXiv:hep-ph/9610451 [hep-ph] . 

[35] E. Witten, Nucl.Phys. B188, 513 (1981). 

[36] I. Balitsky and A. Yung, Nucl.Phys. B274, 475 (1986). 

[37] E. Poppitz, T. Schafer, and M. Unsal, JHEP 1303, 087 (2013), 
arXiv: 1212.1238. 

[38] E. Witten, Nucl.Phys. B202, 253 (1982). 

[39] M. Y. Berry, Proc.Roy.Soc.Lond. A392, 45 (1984). 

[40] M. Stone, Nucl.Phys. B314, 557 (1989). 

[41] F. Haldane, Phys.Rev.Lett. 50, 1153 (1983). 

[42] M. Cristoforetti, F. Di Renzo, and L. Scorzato (Auro- 
raScience Collaboration), Phys.Rev. D86, 074506 (2012), 
arXiv: 1205.3996 [hep-lat]. 

[43] M. Cristoforetti, F. Di Renzo, A. Mukherjee, and L. Scorzato, 
Phys.Rev. D88, 051501 (2013), arXiv: 1303.7204 [hep-lat]. 

[44] H. Fujii, D. Honda, M. Kato, Y. Kikukawa, S. Komatsu, et al ., 
JHEP 1310, 147 (2013), arXiv: 1309.4371 [hep-lat]. 

[45] G. Aarts, L. Bongiovanni, E. Seiler, and D. Sexty, JHEP 1410, 
159 (2014), arXiv: 1407.2090 [hep-lat]. 



